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Abstract. Today, Type la supernovae are essential tools for cosmology, and recognized as major contributors to the chemical 
evolution of galaxies. The construction of detailed supernova progenitor models, however, was so far prevented by various 
physical and numerical difficulties in simulating binary systems with an accreting white dwarf component, e.g., unstable helium 
shell burning which may cause significant expansion and mass loss. Here, we present the first binary evolution calculation which 
models both stellar components and the binary interaction simultaneously, and where the white dwarf mass grows up to the 
Chandrasekhar limit by mass accretion. Our model starts with a 1.6 M Q helium star and a 1.0 M Q CO white dwarf in a 0.124 day 
orbit. Thermally unstable mass transfer starts when the CO core of the helium star reaches 0.53M o , with mass transfer rates of 
1 • ■ • 8 x 10~ 6 M Q /yr. The white dwarf burns the accreted helium steadily until the white dwarf mass has reached ~ 1.3 M Q and 
weak thermal pulses follow until carbon ignites in the center when the white dwarf reaches 1.37 M Q . Although the supernova 
production rate through this channel is not well known, and this channel can not be the only one as its progenitor life time is 
rather short (~ 10 7 - 10 s yr), our results indicate that helium star plus white dwarf systems form a reliable route for producing 
Type la supernovae. 
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1. Introduction 

Type la supernovae (SNe la) are of particular importance in 
astrophysics: They are the major source for iron group ele- 
ments in the universe and are thus an essential contributor to 
the chemical evolution of galaxies (e.g. Renzini 1999 1. And 
their light curve properties allow it to measure their distances 
with an excellent accuracy even out to redshifts beyond z = 
1, which makes them a powerful tool to determine the cos- 
mological parameters (e.g. Hamuy et al. 1996 Branch 1998 
Leibundgut 2000 1. In particular, the recent suggestion of a non- 
zero cosmological constant is in part based on SNe la data 
(Perlmutter et al. [1999 Riess et al. 2000 1. An understanding 
of the progenitors of these supernovae is clearly required as a 
basis for these fundamental astrophysical phenomena. 

However, even though there seems no doubt that SNe la 
are produced by the thermonuclear explosion of a white dwarf, 
it is currently unclear in which kinds of binary systems such 
an event can occur (Livio 2001 1. Among the scenarios which 
have been put forward as possibilities, the so called single 
degenerate scenario is currently favored, where a CO white 
dwarf accretes mass from a non-degenerate companion and 
thereby grows up to the Chandrasekhar mass (e.g. Hillebrandt 
& Niemeyer 2000, Livio 2001 1. However, hydrogen as well as 
helium accretion rates which allow an increase of the CO white 
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dwarf mass due to shell burning are limited to narrow ranges 
(e.g. Nomoto Fuiimoto fT9"S2l Iben & Tutukov 1989 1. 

While hydrogen and helium shell sources are prone to degen- 
eracy effects and related thermonuclear instabilities (e.g., nova 
explosions) at low accretion rates, the helium shell source in 
accreting white dwarf models has been found to be thermally 
unstable even for cases where the electron degeneracy is negli- 
gible (e.g. Cassisi et al. 1998; Langer et al. 2002). This affects 
in particular the potentially most frequent SN la progenitor sys- 
tems, where a hydrogen rich star (main sequence star or red gi- 
ant) is considered as white dwarf companion (e.g. Li & van den 
Heuvei nW71 Hachisu et al.[T333| Langer et al. l2UU0l The hy- 
drogen accretion rates in those systems which may allow steady 
hydrogen shell burning are about a few 10~ 7 M Q /yr. With these 
rates, the subsequent helium shell burning is usually found to 
be unstable (Iben & Tutukov fT^I Cassisi et al. lTWgl Kato & 
Hachisu 1999 1, and no model sequences which cover the major 
part of the white dwarf accretion phase could be constructed so 
far. 

In an effort towards overcoming this shortcoming, we con- 
sider here the evolution of close helium star plus CO white 
dwarf systems. In those, the white dwarf develops only a he- 
lium shell source, thus avoiding complications involved in dou- 
ble shell source models (e.g. Iben & Tutukov I1989I I. Such 
systems form a predicted binary evolution channel (Iben & 
Tutukov 1994 1, which is confirmed — even though for lower 
masses than considered here — by Maxted et al. (2000 1. 
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Further evidence for the existence of close helium star plus CO 
white dwarf systems comes from the recent discovery of a he- 
lium nova (Ashok & Banerjee 2003 Kato & Hachisu 2003 1. 

Simplified binary evolution considerations provided us 
with an estimate for the optimal initial parameters of our model 
system. As a result, we embarked on calculating the detailed 
evolution of a 1.6 M helium star and a 1 M CO white dwarf 
in a 0. 124 d orbit. We introduce our computational method and 
physical assumptions in Sect.|2] In Sect. [3] the evolution of the 
considered binary system is presented. We discuss our results 
in Sect.|H 

2. Numerical method and physical assumptions 

The numerical model has been computed with a binary stel- 
lar evolution code which computes the evolution of two binary 
components, and the evolution of the mass transfer rate and of 
the orbital separation simultaneously through an implicit cou- 
pling scheme (Braun 1997 1. We use Eggleton's ( 1983 ) approx- 
imation for the Roche lobe radius. Mass loss from the Roche 
lobe filling star through the first Lagrangian point is computed 
according to Ritter (119881 . The orbital change due to stellar 
wind and mass transfer is followed according to Podsiadlowski 
et al. ( 1992). The specific angular momentum carried away by 
the stellar wind is computed following Brookshaw & Tavani 
dl993> . Orbital angular momentum loss due to the gravitational 
wave radiation is also taken into account. Opacities are taken 
from Iglesias & Rogers ( 1996 1. Mass loss due to a stellar wind 
is considered as M = 10~ 2 RL/[GM(1 - F)], with F being 
the ratio of photospheric to Eddington luminosity. This mass 
loss rate is based on dimensional arguments and normalized 
for Wolf-Rayet stars (Langer 1989). However, the particular 
choice of the stellar wind mass loss formula does not affect our 
final results significantly, given that the value of the mass loss 
rate is largely determined by F, as the luminosity of the white 
dwarf during the mass accretion phase is close to the Eddington 
limit. Effects of rotation are not considered here (cf. Langer et 
al. 2002, Yoon et al. 2004). For more details about the code, 
see Wellstein & Langer dl999> and Wellstein et al (I200H . 

We start with a zero age main sequence helium star of 1 .6 
M and a 1.0 M white dwarf, in 0.124 day orbit. The metallic- 
ity of the helium star is set to 0.02. In order to avoid the numer- 
ical difficulty in following the initial strong helium shell flash 
in the white dwarf which might be induced soon after the onset 
of mass accretion, the white dwarf is approximated by a point 
mass until the mass transfer rate from the helium star reaches 
about 10~ 6 M /yr (see Fig.[3]in Sect.[3jl. To mimic the heating 
by the initial shell flash, we exchange the point mass by a hot 
and bright (log L S /L Q = 4.175) white dwarf model at this time. 
The initial central temperature and density of the white dwarf 
are T c = 1.8 X 10 8 K and 2 x 10 7 g/cm 3 respectively. Although 
the central temperature may be much lower in reality when the 
mass accretion starts, the initial temperature structure is unim- 
portant for the advanced evolution due to the self-heating of the 
white dwarfs (cf. Sect. |3ji- 

The calculations presented here required high space and 
time resolution: About 500 000 binary models were computed, 
where both stars were resolved into ~ 1000 mass shells. 
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Fig. 1. Evolution of the 1 .6 M helium star in the HR diagram. 
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Fig. 2. Evolution of the CO white dwarf during the mass accre- 
tion phase in the HR diagram. The numbers at the filled circles 
denote the white dwarf mass in the unit of solar mass at the 
given points. 

3. Results 

The evolution of our 1.6 M helium star + 1 M CO white 
dwarf binary system proceeds as follows (see Figs. ^ and [2] 
and Table 1). The helium star undergoes the core helium burn- 
ing for 4.28 million years. After core helium exhaustion, the 
envelope of the helium star expands and starts to fill its Roche 
lobe (Fig.[U when the radius of the helium star reaches 0.61 R . 
At this point, the CO core mass of the helium star reaches 0.53 
M and mass transfer from the helium star commences. Even 
though at this time the helium stars expands on the nuclear 
timescale of helium shell burning, the mass transfer shrinks 
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Table 1. Component masses, orbital period, and central density 
of both stars, for five different times during the mass transfer, 
including the beginning of mass transfer (defines t=0), the time 
of minimum orbital separation (t=73 000yr), and the time of 
central carbon ignition (t=148 000 yr). 
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Fig. 3. Mass transfer rate (dotted line) and white dwarf mass 
accretion rate (solid line) as function of the white dwarf mass, 
for the considered binary model. The upper dashed line de- 
notes the critical mass accretion rate above which the he- 
lium envelope is expected to expands to giant dimensions 
(Nomoto ll982> . For accretion rates below the lower dashed line 
helium shell burning is usually unstable. 

the orbit and thus proceeds on the thermal timescale of the 
helium star. The resulting mass transfer rates are in the range 
1 • • • 8 x 10 6 M Q /yr, as shown by the dotted line in Fig.|3] 

A part of the transferred matter is blown off by the stellar 
wind at the white dwarf surface without being burned into car- 
bon and oxygen, and thus without affecting the thermal evolu- 
tion of the CO core. The rest is transformed into CO in the he- 
lium burning shell source and accreted into the CO core. As the 
mass transfer rate from the helium star increases, the surface lu- 
minosity of the white dwarf reaches about 40% of its Eddington 
luminosity (see Fig.|2li. A radiation driven wind from the white 
dwarf surface is induced, with a peak of 4.8 x 10~ 6 M /yr at 
Mwd- 1-04 M . As the mass transfer rate decreases, the stellar 
wind from the white dwarf becomes weaker. The upper dashed 
lines in Fig. [3] denotes the critical mass accretion rate above 
which the helium envelope is expected to expand to giant di- 




Fig. 4. Evolution of central density and temperature in the 
white dwarf during the mass accretion phase. The dashed line 
gives the locus where the energy generate rate due to the car- 
bon burning equals to the neutrino cooling rate. The numbers 
at the filled circles denote the white dwarf mass in the unit of 
solar mass at the given points. 

mensions, = 7.2 x lO- 6 (M C o/M Q - 0.60) M /yr, based 
on white dwarf models computed with constant mass accretion 
rates (Nomoto ll982> . The helium shell source is expected to 
be stable when 0.4MRh e < M < MRhe- The accretion rate in the 
white dwarf remains in the steady helium shell burning regime 
until M 1 .3 M and the white dwarf mass grows efficiently 
(Fig. [5). Thereafter, the white dwarf undergoes weak thermal 
pulses. 

During the mass accretion, the surface temperature of the 
white dwarf varies from 2.5 x 10 5 K to 12 x 10 5 K and the lumi- 
nosity reaches 2.5 ■ • - 5 x 10 4 L (Fig. [5}, with which the white 
dwarf will appear as a super-soft X-ray source (e.g. Kahabka 
& van den Heuvel 1997 Greiner 2000l. The central tempera- 
ture of the white dwarf decreases initially, mainly due to neu- 
trino emission. When Mwd^I-1 M q , the compressional heating 
begins to dominate, and the central temperature of the white 
dwarf increases. The dashed line in Fig.0]indicates where the 
energy generation rate due to the carbon burning equals to the 
neutrino energy loss rate. The white dwarf reaches this line at 
p c = 1.92 X 10 9 g/cm 3 and T c = 2.71 x 10 8 K, from where the 
central temperature increases rapidly. About 1700 years there- 
after, when T c 2.9 x 10 8 K, the central region becomes con- 
vectively unstable. About 1800 years later, in our last model, 
the convection zone extends to 0.5 M when T c = 4.6 x 10 8 
K(Fig.0. 

The convective URCA process, which is not considered 
in our study, may affect the evolution during this stage 
[1983. 



(Paczynski 119721 Iben 1982 1. It is currently still debated 
whether it leads to heating or to cooling of the convective core 
(e.g. Mochkovitch ll996l Stein et al. I1999i . Thus we have to 
consider the spatial extent and duration of the convective phase, 
as obtained from our models, as uncertain. When the temper- 
ature exceeds 5 • • • 8 x 10 8 K, the nuclear time scale becomes 
comparable to the convective time scale. In this case, convec- 
tion is not able to efficiently carry away the released nuclear en- 
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possible for helium star plus white dwarf systems (Iben & 
Tutukov 1994), this may explain SSSs with various orbital pe- 
riods. For instance, short period systems such as RX J0537.7- 



7034 (3.5h, Greiner et al. 2000) and 1E0035. 4-7230 (4.1h, 
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Fig. 5. Temperature as function of the mass coordinate in the 
last white dwarf model. The dotted lines show the convective 
diffusion coefficient (in the unit of g 2 /s) divided by the square 
of the total mass of the white dwarf. 

ergy, and the mixing length theory becomes inadequate to de- 
scribe the convective energy transport. Therefore, we stop the 
calculation at this point; see, however, Hofiich & Stein (2002) 
for the further evolution. 

In our last model, the helium star has a mass of 1 .08 M Q and 
a CO core mass of 0.64 M . If the white dwarf explodes, the he- 
lium star may survive but suffer stripping off small amounts of 
mass by the supernova ejecta (Marrietta et al. 2000). Its space 
velocity will be close to its final orbital velocity of 330 km/s. It 
will evolve into a massive CO white dwarf of ~ 1 .0 M Q even- 
tually. 

4. Discussion 

The model described above is (to our knowledge) the first 
self-consistent binary evolution calculation which leads to a 
Chandrasekhar-mass white dwarf. Its relevance is in part the 
realistic construction of a supernova progenitor white dwarf 
model, but even more the fact that it provides for the first time 
hard evidence of the functioning of a SN la progenitor channel: 
Our results leave little doubts that some helium star plus white 
dwarf systems will in fact produce a Type la supernova. 

The considered system may have evolved from a wide 8.0 
M giant star +1.0 M G white dwarf system though a common 
envelope phase, which in turn may be the result of two inter- 
mediate mass stars in a close orbit. This indicates a system life 
time of the order of ~ 10 7 — 10 s yr, which is too short to explain 
SNe la in elliptical galaxies. Therefore, our results support the 
idea that at least two SN la progenitor scenarios are realized in 
nature (e.g., della Valle & Livio 1994). Iben & Tutukov ( 1994) 
estimated the potential SN la production rate through white 
dwarf + helium giant binary systems to 1.7xl0~ 3 yr -1 , which 
may constitute a significant fraction of SNe la observed in late 
type galaxies. 

As previously mentioned, SNe la progenitors of the kind 
considered here will appear as a super-soft X-ray source (SSS). 
We note that, since a wide range of orbital separations is 



Schmidtke et al. 1996) can not be easily explained within the 
canonical model which invokes hydrogen rich donor stars (e.g., 
Rappaport et al. 1994 Kahabka and van den Heuvel 1997), ex- 
cept at low metallicity (Langer et al. 2000). Orbital periods of 
the binary system in the present study in the range 2.65 - 2.97 h 
indicate that helium star plus white dwarf systems might be an- 
other natural possibility to explain short period SSSs. 
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